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We present the latest results for the equation of state and the crossover transition in 2+1 flavor 
QCD from the HotQCD Collaboration. Bulk thermodynamic quantities - energy density, pressure, 
entropy density, and the speed of sound - are calculated on lattices with temporal extent N t = 8 in 
the temperature range 140 MeV < T < 540 MeV. We utilize two improved staggered fermion ac- 
tions, asqtad and p4, with the mass for the two degenerate light quarks chosen to be m u( i — 0.1ms, 
corresponding to m K w 220 MeV for the lightest pion. We also calculate observables that are sen- 
sitive to the chiral and deconfing transitions - the light and strange quark number susceptibilities, 
the chiral condensate, and the renormalized Polyakov loop - finding that deconfinement and chiral 
symmetry restoration occur in the same narrow temperature interval. 
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1. Introduction 

Among the most important and fundamental problems in finite-temperature QCD are the cal- 
culation of the bulk properties of hot QCD matter and the characterization of the nature of the 
QCD phase transition. Understanding finite temperature QCD also has direct application in inter- 
preting the results from heavy ion collision experiments. Here, we present a lattice calculation of 
the QCD equation of state, i.e., the pressure, energy density, entropy density, and speed of sound, at 
finite temperature and vanishing chemical potential. We also calculate quantities such as the chiral 
condensate, renormalized Polyakov loop, and the light and strange quark number susceptibilities, 
which are related to the chiral and deconfining aspects of the QCD transition. For a more detailed 
discussion of our results, see [jl|]. 

2. Simulation Details 

This calculation presents results with two different improved staggered fermion actions, p4 
and asqtad, with lattices with temporal extent N t = 8. Both the p4 and asqtad actions eliminate 
the (ff(a 2 ) errors in bulk thermodynamic observables at high temperature, so there are only small 
deviations from the asymptotic, ideal gas limit even at N t = 6 and 8 [0]. Symanzik-improved 
gauge actions are used with both the p4 and asqtad fermion actions. For the p4 action, a tree-level 
improved gauge action is employed. For the asqtad action, a one-loop improved gauge action is 
utilized, with the addition of tadpole improvement in both the gauge and fermion parts of the action. 
Furthermore, to reduce the effects of taste-symmetry breaking, fat-link smearing is implemented. 
The p4 action adds the three-loop staple into the fat link, while the asqtad action adds terms up to 
the seven-link staple which minimize the taste-mixing terms in the fermion action. 

The finite temperature results presented here at N t = 8 are obtained from lattices of size 32 3 x 8 
generated using the RHMC algorithm. These results are presented with previously obtained results 
at N t = 6 [||, f|]. In the case of the asqtad action, some new measurements were also made on 
32 3 x 6 lattices. In addition, "zero temperature" calculations have also been performed on 32 4 
lattices where the gauge coupling j8 = 6/g 2 and bare quark masses are chosen to be the same as a 
corresponding ensemble at finite temperature. 

These calculations were performed with two degenerate light quarks, and a heavier strange 
quark. The strange quark mass is chosen to be close to its physical value, with the light quark mass 
of m„d = 0Am s . The bare masses are fixed so that they lie on a line of constant physics, i.e., the 
hadron masses remain fixed in physical units. In our case, ttik, m Ss , and m n are used to fix the line 
of constant physics. This corresponds to a lightest pion of m % « 220 MeV. The line of constant 
physics differs slightly between the p4 and asqtad actions, with the strange pseudoscalar mass, mj S 
approximately 15% larger in the asqtad case than with the p4 action. In order to set the lattice scale 
in physical units, we use the quantities ro and r\, which are related to the shape of the heavy quark 
potential: 




(2.1) 



We use ro = 0.469(7) fm. in order to translate our lattice data into physical units 
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Figure 1: On the left, the trace anomaly, i.e., (e — 3p)/T 4 for the p4 and asqtad actions at =6 and 8. On 
the right, the same quantity, with the peak of the trace anomaly shown in more detail. 



3. QCD Equation of State 

In determining bulk thermodynamic observables on the lattice, the basic quantity that we cal- 
culate is the trace anomaly &^^(T), 

®^(T) _e-3p_ d / P n 
j4 jA dT \T 4 J ' 

The trace anomaly can be separated into contributions from two pieces. One contribution comes 
from the strange and light quark condensates, and vanishes in the chiral limit, i.e., rhi,th s — > 0, 

-tj^- = -RpR m N? (2m,A (yh^ + m s A , (3.2) 

and a piece that is non- vanishing in the chiral limit: 




Figure 2: On the left, the trace anomaly for the p4 and asqtad actions shown at low temperature, with a 
comparison to the hadron resonance gas. The dashed-dotted/dashed curve is the hadron resonance gas with 
a maximum resonance mass of oimax = 1-5/2.5 GeV. On the right, the trace anomaly at high temperature. 
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Figure 3: Contributions to the trace anomaly from gluonic observables, on the left, and for fermionic ob- 
servables on the right, as defined in the text. Results for both actions at N t = 4,6, and 8 are shown. 



where Rp,R m: and R u are nonperturbative beta functions which describe how the bare parameters 
change with scale: 
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(3.4) 



and the notation A (X) = (X) — (X) T indicates a zero temperature subtraction has been made. 

Figure [l] shows the results for the trace anomaly at N t = 6 and 8. As one can see, there is 
quite good agreement between the two different actions, particularly in the high temperature region 
(T > 300 MeV). Figure |l| also shows the peak region, where the differences between the two actions 
are largest. In particular, the peak for the p4 action is about 15% higher than for the asqtad action. 

Figure ^ shows the trace anomaly at both low and high temperature. Although the p4 and 
asqtad results are generally consistent at low temperature, we find that they (not surprisingly) do 
not agree with those of the hadron resonance gas model, as the hadron spectrum is distorted by 
heavier-than-physical quark masses and taste symmetry violations. At high temperature, we see 
that we have very good agreement between the p4 and asqtad actions, as well as evidence for 
small lattice artifacts. Figure |3| shows the contributions to the trace anomaly from and 0^. 
Although identifying these two terms as "gluonic" and "fermionic" contributions is not accurate (as 
the dynamical quarks mix these two terms), one can see that the gluonic observables dominate the 
trace anomaly. In addition, whereas there is little evidence for cut-off dependence in 0^ for the 
asqtad action, there is significant cutoff dependence for p4. This cutoff dependence can be traced to 
the nonperturbative beta function R m at strong coupling. Indeed, one can see that the cutoff effects 
in 0^ are diminished at higher temperature. 

From the trace anomaly, one can reconstruct the pressure via integration: 

p(T) p(T ) = rT 
ir,, 



f4 



rp4 

y 



J\t>±qw(t>), 



(3.5) 



where To is chosen to be a temperature sufficiently deep in the confined regime so that p(To) is 
small. From the pressure and trace anomaly, it is straightforward to construct the energy density, 
entropy density, and speed of sound. 

Figure || shows the energy density, three times the pressure, and entropy density. Our data 
show that the relative difference between the two actions is no more than about 15% in the low 
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Figure 4: On the left, a comparison of the energy density and three times the pressure for the p4 and asqtad 
actions, atN t =8. The black bar in the upper right hand corner indicates the systematic offset from adding 
the hadron resonance gas value of the pressure at 7b = 100 MeV. On the right, the entropy density for both 
actions at N t = 6 and 8. The crossover region of 185 MeV < T < 195 MeV is also shown on both plots. 



temperature region for T > 150 MeV, and falls to about 5% for T > 200 MeV. For the p4 action, 
the cutoff effects between N t = 6 and 8 are of similar size, while no statistically significant cutoff 
dependence is seen for the asqtad action. 

Figure || shows the ratio of the pressure and energy density, as well as the speed of sound. 
Although there is qualitative agreement between the two actions, the speed of sound significantly 
undershoots the hadron resonance gas value at low temperature. The speed of sound is given by: 

c 2 = dp_ = Jjp/e) 
s de de 



+ ■ 



(3.6) 




Figure 5: Pressure divided by the energy density 
(p/s) and the speed of sound squared (c 2 ) . Data 
points show the ratio of pressure to energy density, 
while colored dashed curves show c 2 obtained from 
the interpolations of s/T 4 and p/T 4 . The dashed- 
dotted line at low temperature shows the result for the 
hadron resonance gas model with m mux = 2.5 GeV 



4. QCD Transition 

As it is suspected that QCD with physical quark masses undergoes a smooth crossover, T c can 
not be sharply defined and may depend on the observable one chooses. In particular, it is unclear 
whether observables related to deconfinement and those related to chiral symmetry restoration give 
similar values for T c . In this section, we examine various quantities related to both deconfinement 
and chiral symmetry restoration. 

One set of observables related to deconfinement are the light and strange quark number sus- 
ceptibilities: 

% q 1 d 2 lnZ 



T 2 VT^dipLq/Tf 



Ls. 



(4.1) 
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Figure 6: On the left, the light quark number susceptibility. On the right, the strange quark number suscep- 
tibility. For both quantities, the Stefan-Boltzmann limit is shown, as well as the crossover r egion (185 MeV 
< T < 195 MeV). 
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Figure 7: On the left, the ratio of the energy density to the quark number susceptibilities. On the right, the 
ratio of the strange to light quark number susceptibilities. The solid curves show this quantity in the hadron 
resonance gas model with m max = 1 .5 GeV (upper branch) and m max — 2.5 GeV (lower branch). 



These quantities give information on the thermal fluctuations of the degrees of freedom that carry 
net quark number. Thus, this observable is sensitive to the liberation of degrees of freedom that 
comes with deconfinement. Figure ^| shows the light and strange quark number susceptibilities. As 
we can see, %i/T 2 rises more sharply than Xs/T 2 , while Xs/T 2 seems to go to zero much more 
quickly. This is explained by the fact that the lightest hadrons that carry light and strange quark 
number are the pions and kaons, respectively. Thus, we would expect Xi/T 2 ~ exp(—m K /T), while 
Xs/T 2 ~ exp(-m K /T). 

In Fig. ^, we see that the light quark number susceptibility tracks the rise in the energy density 
quite well as one moves through the crossover region. On the other hand, e / T 2 x s becomes singular 
in the low temperature regime. This is because the light quark susceptibility, but not the strange 
quark susceptibility, is directly sensitive to the singularity in the partition function in the chiral 
limit. At high temperature, when the temperature is sufficiently high that the quark masses are 
irrelevant, we see that e/T 2 x q agree. The relative suppression of Xs compared to Xi can a l so be 
seen in the ratio Xs/Xu a l so m Fig- where it is consistent with the value obtained with the hadron 
resonance gas model. In Fig. ^ we show the renormalized Polyakov loop and the subtracted chiral 
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Figure 8: On the left, the renormalized Polyakov for the two actions at N t — 6 and 8. On the right, the 
subtracted chiral condensate, A/ s . 

condensate. These quantities are sensitive to the chiral and deconfining transitions, respectively. 
The subtracted chiral condensate is defined as: 

a (T , (M,r-g (W) s ,t 
(VV)/,o ~ % { W),,o 

The temperature band 185 MeV < T < 195 MeV is superimposed on both plots to show that the 
chiral and deconfinement transitions occur in the same approximate temperature regime. 

5. Conclusions 

We have presented here a calculation of the bulk observables (energy density, pressure, entropy 
density, speed of sound) of QCD matter at finite temperature and zero chemical potential. Our 
results indicate that cutoff errors are approximately 15% in the crossover region, and no more 
than 5% for T > 300 MeV. In addition, we have calculated various observables (quark number 
susceptibility, renormalized Polyakov loop, and chiral condensate), which seem to give consistent 
values of T c for the deconfinement and chiral symmetry restoring transition at the lattice spacings 
that we employ. 
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